Relevance of space anisotropy in the critical behavior of m-axial Lifshitz points 
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The critical behavior of d-dimensional systems with n-component order parameter is studied 
at an m-axial Lifshitz point where a wave-vector instability occurs in an m-dimensional subspace 
R™ (m>l). Field theoretic renormalization group techniques are exploited to examine the effects 
of terms in the Hamiltonian that break the rotational symmetry of the Euclidean group E(m). The 
framework for considering general operators of second order in </) and fourth order in the derivatives 
da with respect to the Cartesian coordinates Xa of R™ is presented. For the specific case of systems 
with cubic anisotropy, the effects of having an additional term, X]™=i('^a'/*)^' investigated in an 
e expansion about the upper critical dimension d*(m) = A-\-m/2. Its associated crossover exponent 
is computed to order and found to be positive, so that it is a relevant perturbation on a model 
isotropic in R™. 
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I. INTRODUCTION 

According to the modem theory of critical 
phcnomcna,^'^ systems can be divided into univer- 
sality classes such that the leading singularities of 
thermodynamic quantities of all members of a given 
class are the same. The key concept here is that the 
detailed differences between systems of any given class 
are "irrelevant." The most prominent and best studied 
universality classes are those of the d-dimensional 
n-vector models with short-range (ferromagnetic) inter- 
actions, conveniently represented by the standard (/)^ 
model with Hamiltonian 



n 



i(v0)V^c/>^+ji</.r 



(1) 



Here (p = {(j)a) (with a = l,...,n) is an n-component 
order-parameter field, while f and u are the bare mass 
and coupling, respectively. To access a critical point, f 
must be tuned (corresponding to tuning, say, the temper- 
ature of the system T) to a special value fc so that the 
renormalized r (corresponding to the inverse susceptibil- 
ity of the system) vanishes. The importance of this family 
of models derives from the fact that an enormous variety 
of experimentally studied systems belong to these univer- 
sality classes. Specifically, irrelevant microscopic details 
include the lattice structure, the range of interactions 
(assumed finite and ferromagnetic), or pair interactions 
decaying with a sufficiently large power of the separation 
(e.g., van-der- Waals) 1^ Of course, not all microscopic de- 
tails are irrelevant. Small admixtures of such interactions 
can be treated theoretically as "relevant perturbations" 
to the above universality classes, leading to predictions 
of experimentally measurable behavior of crossovers to 
other classes. 

There are many extensions of the models (QJ, associ- 
ated with systems with more complex types of micro- 



scopic interactions. For example, systems with compet- 
ing interactions may display a richer variety of behav- 
ior. One particular simplified version is the axial next- 
nearest neighbor Ising (ANNNI) modelp*^ in which an 
antiferromagnetic interaction between NNN pairs along 
one of the axes in a simple cubic lattice, in addition to 
the usual nearest-neighbor ferromagnetic interactions, is 
present. By tuning two (or more) control parameters in 
such systems, one can access a Lifshitz point. ^ 

The focus of this paper is the critical properties and 
universality classes of m-axial Lifshitz pointspiSii^^ pos- 
sible in generalizations of the uniaxial (m = 1) ANNNI 
model. To describe these, we split the Euclidean space 
E'* into X with fh = d - m. Let us label the 
coordinates in these subspaces = 1, . . . , m, and 

X/3, P — m + I, . . . , d, respectively, and introduce the no- 
tations da = d/dxa and dp = d/dxp. Then the Hamil- 
tonian of these extended models reads 

pom 

Hiso = Jd'x (^"'^)' + 2 ^ (^'3'^)' 



+ 



CTi 



(2) 



Provided the microscopic aspects (e.g., d > d^{m, n), the 
lower critical dimension hera^) allow the system to be 
tuned to Lifshitz points, they occur at critical values plp 
and tlp of p and f. Analogous to the above case, both 
of the renormalized parameters vanish at these points: 
p and T ~ Q. For m = 0, the model |2Jl reduces 
to the standard isotropic cj)'^ theory of Eq. with no 
Lifshitz points. At the other extreme, m = d, the system 
displays an isotropic Lifshitz point. 

Although the model ||2J) was introduced more than 25 
years ago, its systematic investigation beyond Landau 
theory via modern methods of field-theoretic renormal- 
ization group has just begun liSiiiiiSiiiiiii^ Within the 



2 



framework of an e expansion (e = d*{'m) — d, d*{m) is 
the upper critical dimension 4 + m/2, and < m < 8), 
early studies were either restricted to special values of 
m and a subset of critical exponentaSiiSiii or else pro- 
duced resultSiiZiil to O(e^) in conflict with those of Sak 
and Gres^ and more recent workiifl*ii*iSiiiiiiii Only 
recently has it become possible to master the enormous 
technical difficulties one encounters beyond the one-loop 
approximation. The full two-loop renormalization group 
(RG) analysis yielded results for all exponents (critical, 
crossover, and correction-to-scaling) to order for all 
values of 77j,i2iiiiiii^ An alternative picture of the Lif- 
shitz point has been advocated by Leiteii^ This has been 
critically assessed in Ref. Il^ 

Let us also mention some earlier work on modifica- 
tions of the model ||2J). Hornreich^^ investigated the ef- 
fects of contributions breaking the 0{n) invariance of the 
Hamiltonian, using a one-loop approximation. Folk and 
Moser— studied Lifshitz points in systems with short- 
range and uniaxial dipolar interactions such as uniaxial 
ferroelectrics. 

The purpose of the present paper is to examine the le- 
gitimacy of taking the fourth-order derivative terms of 
the Hamiltonian (0) as isotropic in the subspace R'". 
Made essentially for the sake of simplifying the com- 
putations, this assumption of "m-isotropy" is question- 
able, since the discrete lattice symmetries at microscopic 
scales are unlikely to respect full rotational invariance. 
Of course, we must account for these underlying symme- 
tries at the continuum level when appropriate Hamilto- 
nians are considered. Now, in the long wavelength limit, 
isotropy can be restored by appropriate rescaling of the 
axes at the level of second order derivatives. However, 
there is no such luxury in general, at the higher orders. 
Hence, the replacement 



CTl 



a=l 



should be made in Eq. where T is a hnear combina- 
tion 



I aia2a3CK4 



(4) 



of tensors T''-* compatible with the symmetry of the mi- 
croscopic model considered. Here the summation conven- 
tion is used: The doubly occurring index i as well as all 
pairs of a indices are to be summed over. In general (the 
"m-clinic" case, a generalization of the familiar triclinic 
one for m = 3), there are 



m + 3 
4 



(5) 



such tensorsi22i24 Instead of dealing with the Um coeffi- 
cients (7i, a convenient set are the rim — 1 dimensionless 
coupling constants 



along with ai . 

In this paper, we will focus our attention on a familiar 
example: The symmetry associated with an m-cube, i.e., 
cubic anisotropy. Besides the totally symmetric tensor 

(7) 

we have only another one, namely, the cubic 

^aiQ2fi3CK4 ~ ^aia2 ^a2fi3 ^(^304 ■ (8) 

Thus Eq. 10} reduces to 

and the Hamiltonian becomes 



^2 / ,d 



(10) 



Note that this model should represent the universality 
class of a simple generalization of the ANNNI model, i.e., 
from the uniaxial Ising (m — n = 1) to the m-axial 0{n) 
case. Specifically, consider a simple cubic lattice Z'* with 
classical n- vector spins Si of unit length on its sites i. 
Assume that the spins are coupled in an 0(n) symmet- 
ric fashion, but with different characteristics within the 
two subspaces R™ and R™. In the former, suppose the 
interactions are like those in the ANNNI model: nearest- 
neighbor ferromagnetic (of strength Ji > 0) but sec- 
ond neighbor antiferromagnetic (strength J2 > 0) along 
each of the m principal lattice directions. In the comple- 
mentary subspace, let the interactions be only nearest- 
neighbor ferromagnetic (of strength J3 > 0). The lattice 
Hamiltonian Tiiat is explicitly 



kBTHla 



E 



(11) 



Denoting the Fourier transform of Si by Sq, we recast 
this expression in Fourier space: (qa, <![})■ The first line 
yields a contribution cx X]™=i [•^a cos(2(7q) — JiCos{qa)] 
to the coefficient of |sqp. The Lifshitz point can be 
accessed, at this naive level, by tuning the 0{q^) term 
to vanish, i.e., Ji = 4J2. Meanwhile, the 0{q^) term 
is precisely of the form of the a2 term in the coarse- 
grained Hamiltonian (|10|l . Though this procedure does 
not directly yield a cri term, there are two good reasons 
that such a term is unavoidable. Firstly, we generalized 
the ANNNI model in the simplest possible manner: All 
further-neighbor interactions couple only spins along the 
principal directionSfS^^SS, Had we introduced NNN bonds 
along diagonals (i.e., i — j — ±6^^ ie^^), there would be 

|s„P, which in- 



a contribution of the form / 



(6) volves both the symmetric tensor ||7J) as well as the cubic 
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one (|SJ. Secondly, the isotropic coupling will be auto- 
matically generated when short wave-length degrees of 
freedom are integrated out, as will be shown in the RG 
analysis below. Thus, we expect that a wide class of lat- 
tice models similar to will fall into the universality 
class described by the Hamiltonian (|10|l . 

Unless stated otherwise, we will restrict our attention, 
for simplicity, to this case (i.e., ©, tensor with cubic 
symmetry) and study only the Hamiltonian H1U|I . Our 
goal is to examine the effects of this type of anisotropy 
on the isotropic m-axial Lifshitz point. Generalizing a 
two-loop RG analysis of the latter casefiSiii we will show 
that the cubic anisotropy cx (T2 is a relevant perturbation, 
at order e^. 

In the next section we present the formal framework 
for renormalization of the m-anisotropic model with gen- 
eral tensors of the form , including the associated RG 
equations. In Sec. Ill Bl we specialize to the case (fTU)) with 
only a cubic anisotropy. Since the anisotropy of interest 
appears in the momenta of a two point vertex function, 
a two-loop computation is necessary. More explicit re- 
sults, to first order in the cubic anisotropy are provided, 
so that its effects on the RG flow near the isotropic fixed 
point, as well as scaling properties, can be investigated. 
For general values of m and n, the e expansion, to O(e^), 
of the associated crossover exponent, if2{n,m,d), is ob- 
tained in terms of integrals over a single variable. In 
Sec. mil we compute these integrals, analytically for the 
special cases of m = 2,6 and numerically for a range of 
other m's. An estimate of 1^2(1,2,8) is presented. Con- 
cluding remarks are reserved for Sec. IIVI Finally, there 
are three appendixes to which some details of our calcu- 
lations have been relegated. 



II. RENORMALIZATION GROUP ANALYSIS 

A. General anisotropy: Renormalization and RG 
equations 

To renormalize our theory with general "m- 
anisotropy" we straightforwardly extend the 

considerations for the m-isotropic model Q in Ref. 0. 
For the details in the analysis, we will follow the 
conventions and notations of Ref. Il3l Here, we have 
n„i variables, ai, whose scaling dimensions vanish at 
the Gaussian fixed point p = t = u = Q. There are 
two consequences: Associated with each of the cji is a 
renormalization factor Z„^. Further, these quantities, 
as well as all other renormalization factors, become 
functions of dimensionless coupling constants, 

namely, the renormalized four-point coupling u and 
nm — 1 renormalized counterparts of the bare variables 
Wi [Eq. ©]: 

Wi = (Ji/ai, i^2,...,nm ■ (12) 



Accordingly, we reparametrize the theory as 





= Z^. (U, W) (Tj , 


(13) 






(14) 


/O ON 0-1/2 

(P-PLP) cr^ 


= p,Zp{u,w)p , 


(15) 


T - T-LP 


= p^ Zt-{u,w) [t + Ar{u,w] 


ip2],(i6) 




— Zu{u, w) u , 


(17) 



where is a momentum scale, Fm,iL denotes the normal- 
ization factor 

p _ r(l + e/2)r^(l-e/2)r(m/4) 

(4^)(8+m-2e)/4r(2-e)r(m/2) ' ^ ' 

and w stands for the set {wi. ■ ■ ■ , Wn^} of — 1 vari- 
ables. Following Ref. 1271 we have included a renor- 
malization function Ar{u,w) to absorb momentum- 
independent poles proportional of the two-point vertex 
function. 

The fact that the theory must reduce for w = Q to the 
m-isotropic one implies the relations 

Z,{u,w = Q) = Zf^{u) , L = (t),u,T,p, (19) 
Z„,{u,w = Q) = Zf'iu) , (20) 

where the Z factors marked by the superscript "SD" are 
those of Ref. [13 The function Ar (u, w — 0) has been 
computed to one-loop order in Ref. f27l Its explicit form 
will not be needed in the sequel. 

Turning to the RG equations, we use the notation 9^|o 
for p derivatives at fixed bare variables (u, ai, f and p), 
and define the (3 and exponent functions 



(3i^ = pdf^li^ K , K = u,T,p,a^ , (21) 

rix{u,w) = pd^l^lnZx , A = (/), m, t, p, CTj . (22) 

The functions rjx depend only on u and w. Since we use 
minimal subtraction of poles, they are even independent 
of e. In terms of these variables andSi 

br{u,w) = Arlpd^l^lnAr +f]T -iTJp] , (23) 

the P functions can be written as 

(3u = -[e + ri^iu,w)]u , (24) 

(3r = -[2 + T]r{u,w)]T - p^ br{u,w) , (25) 

I3p = -[l + Vpiu,w)]p, (26) 

Pa, = -?7<T. {u, W) CTj . (27) 



Also, since we fixed all renormalization factors Zx such 
that the regular part of their Laurent series in e is ex- 
actly unity, the associated rjx functions are related to the 
residues of the Zx 's via 

ijx^ -uduResZx . (28) 
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In terms of the operator 



(29) 



the RG equations of the renormaUzed A''-point cumulant 
functions Gr^n{x) = ( IljLi </'aj,ren(a3j))''"™ and the cor- 
responding vertex functions Fr^-* read, respectively, 



N ■ 



g(^) = 



N 



r(w) = 



(30) 

Being dimensionless, w will appear in the Z's to arbi- 
trary orders in general, even though we are dealing with 
the systematics of an expansion in powers of u (the loop 
expansion). However, our principal goal here is a local 
stability analysis of the model HlOfl about the isotropic 
fixed point 



0) 



(31) 



where u* is the nontrivial zero of /3„(e,u,K; = 0) the 
exphcit form of which, up to O(e^), is given in Eq. (60) 
of Ref. To this end, we can linearize about Vig^- 
Hence it will be sufficient to compute the counterterms 
to first order in w. 



B. Cubic anisotropy: RG flow and scaling 

Given the general framework above, we turn to ex- 
plicit results for the particular case with only a cubic 
anisotropy. With just cti and a2, we have only one w 
(and no need for the set w). Further, we need only terms 
linear in w. Referring the reader to Appendix 1X1 for the 
computational details, we note here that the pole terms 
of the two-loop graph ^ ^ yield the renormalization 
factors 



1 



3 12 8 - m) 



(32) 



j^(m) — 36i^{m)w — + 0{w ,u ) , 



3 96 m(m-)- 2) 



(33) 



jcr(TO) — 36 Vi (m) w hO(u;^,u'^) 

e 



and 



ZZ = 1- — 4v,(m) v?_ 

'^^ 3 TO(m -f- 2)(m -|- 4)(m -|- 6) e 

2 „,3\ 



0(wu , u 



(34) 



Here, j^{m) and 7'^ (to) are single- variable integrals en- 
countered in Ref. \l3L their definitions noted for the read- 
ers' convenience in Appendix m Eqs. (|X28|) and l|J29|l . 



Though more complicated, i^{m), (jn), and (to) are 
analogous integrals, defined in Eqs. (|A24|I - HA26|I . 

From the structure of these Z factors, it is clear that, 
at this two-loop order, both rj^-^ and 770-2 are proportional 
to ^^^M^ with coefficients linear in w. As a result of 
Eq. (|27|) . and keeping only terms to first order in (72, we 
find the associated (3 functions to be of the form 



/3(T2 

where 



= -2- 





K12 \ 




\ 


K22 / 





[l + 0{w)] . 
(35) 



_ j4,{m) 



K12 = - 
K22 



12 (8 -to) 96TO(TO-f2) 



8 - TO 8m(TO + 2) ' 



12 I 



to) m(TO 2)(to + 4)(m + 6) ' 



are constants, independent of the couplings. 

Since our main interest is the neighborhood of V^^^^, we 
need to evaluate u in this equation only at the TO-isotropic 
fixed point 



6e 

n + 8 



(36) 



Given the form of the matrix in Eq. (|35|l . the eigenval- 
ues are trivially obtained, the first of which is just 77* in 
Refs. Eland El 

24(n-h2) 2 r^r 3^ 

- ^^ = -7;rT8F^"^' + o(e3), 

, _ 24^n + 2) 2^nf3^ 

Associated with these are, respectively, the (linear) scal- 
ing fields: 



(Ji -\-h1J2 and (72 



(37) 



where h = Ku/Kn — K22- Near 7^;*^, we may drop the 
irrelevant contributions proportional to m — it*, so that 
the flow equations 



(38) 



are solved by (7i(^). Imposing the initial conditions 
CTi(l) = (7i, these take the asymptotic forms 



« (cTi +6(72)^""^ -6CT2r''- 
a2(£) « (72 r"^^ . 



(39) 
(40) 



Thus, to the order of interest, the anisotropy w leads to 
a dependence on the following ratio of running variables: 

a2(^)M(^) « ((72/(71) rf''^^-'';) (41) 
= [w + 0{w'')\ £^("^2-";) . (42) 
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As we will show, rj*^ — rj* is positive, so that the effect 
of w is more significant in the infrared limit {£ — *■ 0). 
Therefore, we introduce the anisotropy crossover expo- 
nent 



(43) 



which governs the scaling behavior of w with t: w ^ 
T~'^^ . Since the rj's are already of O(e^), we may insert 
the zeroth-order value for 1^12 (i.e., 1/2) to obtain 



1 



'P2 



{n + 8)2 m{m + 2) 
48 10-2(771) 
(777 + A){m + 6) 



jo (777) 



8 



(44) 



As a consequence of the contribution proportional to 
of [see Eq. (|25|) ]. the variable t is not a scaling 
field. Proceeding similarly as in Ref. i27t we can define a 
nonlinear scaling ficld^— 



(45) 



with the asymptotic scale dependence gri^) g^^ 
where the ellipsis stands for terms of higher order in w. 

Utilizing the above results, one can generalize the con- 
siderations of Refs . IT^ andl2^ in a straightforward fashion 
to obtain the scaling forms of the renormalized iV-point 
cumulants G^^^: 



G^^^ [{xq, x^}; gr, p, (Ti,cr2,u, p.)] 
9t 



f y/pXg pXf3 \ p_ W 



5t J 5t gf 

where Ag is the scaling dimension of G^^^: 

Ag = (iV/2) [d-2 + rji2+m{e- 1)] 



(46) 



(47) 



With the exception of ip2, explicit expressions for all ex- 
ponents, up to O(e^), are given in Refs. IT^ and IT^ In 
the next section, we present a brief summary of the steps 
for computing ip2. 



III. EXPLICIT RESULTS FOR THE 
ANISOTROPY CROSSOVER EXPONENT 

From Eq. H44I) . we see that only two integrals are 
needed for finding the exponent ip2, namely, jaifn) and 
7ct2 (w). As shown in Appendix IXI all integrations can be 
performed analytically, except the one over the scaling 
variable 



ixpxfj) 



1/4 



(48) 



for which we resort to numerical means. The first, jo, is 
familiar from Ref. fisl recalled here in Table 



The second involves, in general, a product of three hy- 
pergeometric functions. As outlined in Appendix E of 
Ref. 0, we write 



7^2(777) 



dvlrn{y) 



dv Irniv) + Rrniva) , (49) 



splitting the integral into a contribution from a finite 
interval (0,uo) and a remainder Rm{vQ). The point vq 
was chosen so that standard numerical integration rou- 
tines (specifically. Mathematical^) yield sufficiently ac- 
curate results for the first term and a few terms of the 
asymptotic expansion for suffice for evaluating Rm 
approximately. In practice, we chose vg ~ 9.5 and just 
the leading term of the asymptotic expansion. The latter 
can be computed analytically and leads to: 



Rmivo) 



3V7r 2«-"(777 - 2)2 

r(f)r(i4 



r 2 



(50) 



Combining this with the numerical integration over 
(O,wo)j we arrive at the values of v, (777) displayed in 
Table |I[30 

As discussed in Ref. 0, the cases m — 2 and ?ti = 6 are 
quite special. The scaling functions from which Im{v) is 
formed — and hence Im {v) itself — reduce to elementary 
fimctions. (From another perspective, their asymptotic 
expansions terminate at low ordersiiSii^ See Eqs. ljCl|l . 
(|C7|I and (|B8|I . (|B9|I . Indeed, the approximation l(5?Hl 
even vanishes for 777 = 2.) As a result, i^^ can be com- 
puted analytically. The results (see Appendix EJ, 



2 / 4 59 

V.(2) = — , 7o2(6) = 20 Sin-- — 



27 



27 



(51) 



are in conformity with the values quoted in Tableland 
provide useful checks of our numerical procedure. For 
completeness, let us also recall the analytic values 



128 448 
J<t(2) = — , = 



(52) 



from Ref. [H 

Inserting the values from TableUinto Eq. H44|) , we write 
the anisotropy crossover exponent as 



^2 



27(77 + 2) 
(77 + 8)2 



G(777) e^ + 0{e^) 



(53) 



Here the C{m) are the coefficients of the crossover 
exponent ip2 for 77 = 1 . They are listed in the last column. 

For the special cases where analytic results are avail- 
able, we have 



G(2) 



1 

324 



and 



C(6) ^ - 



4 4 19' 

3 ^^3^81 



(54) 



(55) 
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As we see, in all cases, 1^2 is positive at this order. Thus, 
we conclude that the isotropic fixed point V*^^ is unstable 
against perturbations from cubic anisotropy of the form 
included in the Hamiltonian pU|l . 

TABLE I: Numerical values of the integrals ja{rn), ia^{m), 
and the coefficients C(m) of ip2 introduced in Eq. 1531 . 



m 






C(m) 


2 


4.74074 


0.074074 


0.00309 


3 


10.804 


0.24682 


0.00380 


4 


20.067 


0.6175 


0.00444 


5 


32.95 


1.279 


0.00501 


6 


49.7778 


2.325428 


0.00552 



Of specific interest is the scalar biaxial case m = 2, 
n = 1, corresponding to a biaxial generalization of the 
three-dimensional ANNNI model discussed in the Intro- 
duction. The upper critical dimension being 5, let us set 
e = 2 to consider a physical system in d = 3: 

(^2(n=l, m=2, d=3) 0.0123 . (56) 

81 

Since this exponent is so small, we see that, unless the 
anisotropic amplitudes are large, we must be extremely 
close to the Lifshitz point (by careful tuning of the two 
control parameters) in order to detect any serious de- 
viations from the critical behavior in the class of the 
isotropic fixed point V'^^^. 



form included in the Hamiltonian (|10|l will be generically 
present even in systems with lower symmetries, leading to 
crossover in general. Since all previous investigations^ 
of rn-axial Lifshitz points (with m > 1) are based on 
m- isotropic Hamiltonians, our conclusion is that, un- 
less the microscopies of a system enforces rotational in- 
variance in R™, the critical properties will not fall into 
the universality classes found so far. Nevertheless, we 
should caution that, though the crossover exponent is 
positive, its numerical values are relatively small. For 
example, in the case most likely to be physically acces- 
sible (d = 3, TO = 2, n = 1), this exponent is only 1/81. 
Though measurable deviations from the isotropic class 
may be difficult to detect for real systems, it should be 
interesting to test our predictions in Monte Carlo simu- 
lations of suitably designed lattice models. 

A more interesting question is, given the RG flow is 
away from the TO-isotropic fixed point, whether there is 
a new stable fixed point or not? Of course, if the former 
is true, then we have a new universality class. However, 
preliminary studies indicate that the (3 function for w has 
no zero, for any finite, non-vanishing w. Worse, the flow 
seems to run indefinitely towards Hamiltonians associ- 
ated with singular propagators (i.e., parts quadratic in 
</> being no longer positive definite), so that 0{q^) terms 
will be needed to stabilize the theory. In addition, we 
might speculate that the system would undergo a first 
order transition, albeit a weak one. Clearly, further in- 
vestigations are necessary before definitive conclusions 
can be reached. 



IV. CONCLUDING REMARKS 

We studied the critical properties of m-axial Lifshitz 
points in systems with spatial anisotropic interactions. 
Specially, for an TO-dimensional subspace of in which 
a wave- vector instability occurs, we considered the effects 
of arbitrary fourth-order couplings of the form 

'^ia2Ct3Ct4 QaiQa2^0(3Qa4 ' t (^'^) 

with 4> being an n-component order parameter. Unlike 
previous studies of "anisotropic" Lifshitz pointSfSi we are 
not concerned with couplings that break the 0{n) sym- 
metry of the order parameter. In this sense, our (f) is 
more appropriate for, say, multi-component alloys than 
for spin systems. Especially for the n = \ (Ising) case, we 
showed how anisotropy of the form l|57|l naturally arises 
from a generalized ANNNI model. In the case we explic- 
itly considered only a cubic anisotropy, corresponding to 
a tensor (|SJ), was present [cf. Eqs. © and HlU|l ]. 

Using field-theoretic renormalization group techniques 
and expanding about the upper critical dimension, we 
found that this is a relevant perturbation for an isotropic 
TO-axial Lifshitz point. To O(e^), the crossover expo- 
nent has been computed for a range of m, with analytic 
forms for the special cases of to = 2, 6. Though we have 
not obtained similar results for general anisotropic inter- 
actions, there is no doubt that cubic anisotropy of the 
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APPENDIX A: CALCULATION OF THE 
TWO-LOOP GRAPH OF T^^^ TO ORDER era 

Let us denote momenta in M'' by q = {k,p) with com- 
ponents k e K™ and p £ M™ {fh = d — m). As we will 
be dealing with cubic spatial anisotropy extensively, let 
us define 

m 

fc'^E^- (Ai) 

a=l 

which is to be distinguished from the usual fc^ = {k-kY . 
At the Lifshitz point t = p = 0, the Fourier transform of 



7 



G{x), the free propagator in position space, is 



G(q) = G(fc,p) 



(A2) 



To two loops, r(^)(<7), the Fourier transform of the bare 
vertex function F*^^^ , is given by 



G{q) 



^^GHq) + Oiu'), (A3) 



where G'^(g) denotes the Fourier transform of G^{x). In 
addition, at this order of the e expansion, we will need 
only the residues of the simple pole of F'^' at e = 0. 
Thus, all coefficients of may be evaluated at e = 0, 
or d = d* = 4 + m/2. 

Let us compute the pole terms of G^ to first order in 
w = (72 /(Ti- Since the dependence on cti follows from 
dimensional considerations, we temporarily set ai = I. 
Expanding in w, we have 



G3 = 



G3 



9,;,G3 



+ 0(u'') 



(A4) 



The pole contribution of the first term on the right-hand 
side has been computed in Ref. ^3- Expressed in terms of 
the integrals j^{m) and (m) of Ref.[l3 [see its equations 
(43), (44), and (46)], the result is 



G3 



w=0 



'2(8-m)^ 



16 r7i(m + 2) 



+ 0(6") 



(A5) 



To compute the w derivative appearing in Eq. (|A4p . we 
start from its position-space representation, G^{x), and 
obtain 



[d^^G\x)]^ 



-iGl{x)Y,di{G^*G^){x) . (A6) 



Here Gq{x) = G{x)\^=q, and (Gq * Gq){x) means a 
convolution in position space (i.e., Fourier transform of 

[G,{q)?). 

Next, let us exploit the scaling forms of both Go and 
Go * Go. Defining the radii 



ex. 1 



the scaling variables 

Va = Xa R~^^^ , Z = V^Va = , (A8) 

and the vectors v = {va} and e = {xp^/R^ we write 



Go(a;) = $™,d(t;) 



and 



(Go*Go)(a;) = i?^r,„^<j(z) 
Here, the scaling functions are 



and 



where /^"^ = (27r)-™/d"A: and = (27r)-"/d"V. 
Their explicit expressions in terms of Taylor series and 
hypergeometric functions, as well as the asymptotic ex- 
pansions for large z are given in Appendix IBl 

Inserting the above scaling forms into Eq. (|A6p and 
performing the Fourier transformation, we encounter 
integrations over x. Using hyperspherical coordinates 
{vVR, f2„j) in R™ and (i?, 17^) in E™, let us denote an- 
gular averages by 

j"" = Sa'Jf{nD)dnD (A13) 

where Sd = f dflo = 2tt^/^ /r{D/2) is the surface area 
of a _D-dimensional unit sphere. In the radial integra- 
tion dR, the distribution R~^~^^'^ is found to appear. 
Employing its Laurent expansion^S 



(m) pirn) ^ik-v+ip-e 



(m) j-{m) ^ikv+ipe 



(fc4 + p^y 



(A9) 



(AlO) 



(All) 



(A12) 



R^^+^' = — S"{R) + 0(e°) , (A14) 



(A7) we find 



a,;,G3 



4e 



+ 0(6°) 



R=0 



The right-most part of Eq. (jAlSp can be rewritten as 

■m 

a—1 a 

where y(*^)(z) = d-^Y (z) / dz'^ . 



(A15) 



(A16) 
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Note that the function of R that must be differentiated 
in Eq. (| A 15(1 has an expansion in even powers of ^/R be- 
cause the angular integration in J d"^v yields zero for 
the coefficients of all odd powers. Thus only the contri- 
butions oc produced by the two exponential functions 
contribute to the pole term on the right-hand side. 

For the angular averages, we may exploit the charac- 
teristic function 



ik-v 



^4^£!r(£+f) 



(A17) 



to find 



ing (2A'^ — 1)!! — 1 pairings of the indices. This can also 
be thought of as 2A^-point function of a Gaussian the- 
ory whose propagator between two "points" and aj is 
given by Sa^a, ■ 

Utilizing these results, we easily find 



dR^ 



R=0 



in 



E—rm 3 
< = ——7; : 
m + 2 

a 



(A19) 



(A20) 



■m i,2N 



^a.ia.2 ' ' ' ^Ct2N - lO'-2N 



(A18) 



m(m + 2){m -f 4) • • • [m 2(iV - 1)] 
In the numerator the ellipsis . . . stands for the remain- 



and 



_R=o 4:m{m + 2) 



(A21) 



92 



a 

Combined with Eq. (|A15|I . these results yield 



3(to + 8)k^ + 8 
_R=o 4m(m + 2)(m -|- 4)(m -|- 6) 



m.e 

w=0 e 



18 i^{m) 2 9Vi(m) 4 



■ p 

8 — m 4r7i(m -|- 2) m(TO 4- 2)(to + 4)(m + 6) 

where we have introduced the integrals 



24v,(m) 



+ O(e0) 



(A22) 



(A23) 



Jo 



m 



4 v{™) I 1\ 



m + 2 "'''^^ 



(A24) 



and 



iai (fn) = B„i / dv V 
Jo 



m+3 



mYZ (v^) + 4 Y-,, (v^) + / + Y^. 
rn,d \ / m.d \ ; (m -|- 4) (m -t- 6) 



(A25) 



z., (m) ^ / dv <^l,,, (v) Y^, (v^) 



(A26) 



with 



r(2-TO/4)r2(m/4) 



f2 

^m,0 



(A27) 



For completeness, we recall the definitions of 



and 



Jo 



(A28) 



Jo 



(A29) 
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APPENDIX B: THE SCALING FUNCTIONS $™,d(w) AND Y^,d{z) 



First, let us recall the properties of the scaling function $m,ti(w), which were established in Refs. andll^ Its 
Taylor expansion is respectively, 



i=0 ^--^12 ' 2 



4 y 



T 



(Bl) 



It is possible to express the power series in closed form by exploiting a relation of the generalized hypergeometric 
functions 1F2 (by summing the contributions with even and odd values of i separately), namely. 



- 1 r(a + |) 



T{a) 



1 y2 

1^2! a; 



r(a + i) 



The large 77 properties of 1F2 lead us to the asymptotic expansion 

21-" TO -2 



r(5 + i) 



1 3 

2' 2' 



ly 
2' 4 



$m,d*(w) 



.00 7r(6+™)/4 r(i + f ) 



i + o(u-^) 



(B2) 



(B3) 



which will be used in numerical evaluations of i^^ (to) described in Sec. IIIII 

The scaling fmiction YmA'^^) "^^^ ^-l^o studied in Ref.0: It appeared in the calculation of another scaling function, 
Q{v), that was needed for the coupling constant renormalization at the two- loop level. Its Taylor expansion, from 
Eq. (D4) in this reference, reads 



■(-i)' 



As for <& above, we can write a closed form for this function in terms of two hypergeometric functions: 



r(i-rf)/2 

24+m 



1F2 



el m 



2' 2 4 '647 4r(l + ^ 



^1. , m\ 1^2 t: - 77, 7:, i 



2' 2' 



m z 
T' 64 



(B4) 



(B5) 



Let us note that, though the series HB4() starts with a term which has a pole in e, all its derivatives are regular at 
e = 0. In particular, for the calculation of the integral ic2{ni), we need only Y^^J,[z) at the upper critical dimension 
d* = 4 + to/2. Thus, 



2l2+m_(6+m)/4 p\ y 

1=0 ■ V2 ' 4 



^i r(2 + f) (-1)^ 



1 



2l4+m 7^(6+m)/4 



r(i + f) 



i^^2 2;-, 



15 m 



- 



2' 2 4 '64/ r(3 + 2i 



1F2 



5 3 m z 
2' 2'^^T'64 



In the special cases to, = 2 and to = 6, the hypergeometric functions here reduce to the elementary functions: 



1 - e 



-z/4 



32 384 



and 



(47r)2z4 

48 (-320 + z2) + [15360 + z (I6 + z) (240 + 12 z + z'^)] 

(4 7r)3 64z6 

I 



(B6) 
(B7) 

(B8) 
(B9) 



For the integral ^0-2 (to), we will need the asymptotic expansion of y^^^'^l(z) 



(BIO) 



22+m 7^(6+m)/4 1 



l + 0(z-2) 
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Let us however note that, for arbitrary values of to, 
derivatives of Y^.d* (z) can be expressed completely in 
terms of ^m,d'{v) (with v = Y^)fi2ii^ since their Taylor 
series expansions differ only slightly. In particular, using 



we get the general relation 



e 

'--2 



2 2/ V2 2 



(B12) 



Recalling that Ym^d has a pole at e = 0, we verify that 
Ym,d{z) = -$m,d.'(0)/(2e) + O(e0). Thus, in the limit 
e — > 0, we find an exceedingly simple relationship be- 
tween Y^ ^ and ^jn.d- 



Ym.d* (^) 



1 r 

2i 



(B13) 



In other words, the scaling function '^m,d* completely 
determines all derivatives of Ym^d' , a fact we will exploit 
in the following Appendix. 



APPENDIX C: ANALYTIC RESULTS FOR m = 2 

AND m = 6 

The special case to = 2 is the simplest one since the 
scaling function <I> acquires an extremely simple form at 
the upper critical dimension: 



1 



(47r)2 



(CI) 



For the function 51 using the relation (|Ijl3() . we thus 



get 



1 1 



(4^) 



2z 



-z/4 _ 



(C2) 



Even more simplifications of integrals are realised if we 
use the representation 



C dte-^'/^^ -(1 
Jo z 



(C3) 



so that 



.(if)/ N 1 2 



(4^)2 45 7o 



(C4) 



(Bll) i^^^ we find 



Inserting this expression into Eq. (|A26|I for the integral 



/>00 /"l 

V,(2) = 2 dCCe-^'^ / dtt 
Ja Jo 



(C5) 



To arrive here, we changed the integration variable to 
^ = z/4: for convenience. Carrying out the trivial inte- 
grations, we get 

i.,{2) ^2-41 J' dtt^t + 2)-'^^. (C6) 

The special case to = 6 is, as usual, somewhat more 
involved. This is due to a more complicated functional 
form of the corresponding scaling function Using ^ = 
z/4 once more, we have 



$6,7(V4C) 



1 1 r 



2(47r)3 



1-(1 + C)e-'^ 



(C7) 



Again, we exploit an integral representation like Eq. (jC3|) 
— now, for both functions $6.7 and 7. Thus, (|C7|I can 
be written as 



$6,7 (V4C) 



1 



dt t e"*^* . 



2(4^)3 7o 

Inserting this into the relation ljB13p . we obtain 
1 



(C8) 



and 



1 i-i 

2 / J -zty/i 



16 (477)3 



dtt dye 



1 /.I 

5/ ^„„3„-zty/4 



dtt'' dyy-^e- 



4^4^)3 7o J, 
The desired integral (6) thus becomes 



(C9) 



(CIO) 



i^^(6) =4 dCC^ daae-°"^ dppe''^'^ / dttN dyy'^ e'^'y 
Jo Jo Jo Jo Jo 

1 



4-6! / daa [ dj3 fi [ dtt^ [ dyy^ 
Jo Jo Jo Jo 



/ 4 59 
ia + P + tyy-H''^^3-27 



(Cll) 



The integration over C was trivial, and the last result was of MathematicAiSS 
obtained by simple repeated integrations with the help 
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